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Abstract 

A new neutron lifetime experiment employing frozen Fomblin has produced a result in significant 
disagreement with previous experiments that used liquid Fomblin near room temperature. This 
new experiment is subject to very few corrections, so the source of the discrepancy remains to 
be identified . Here we theoretically investigate several possible systematic effects for near-room- 
temperature storage experiments. By considering the combined effect of gravity with the usual 
ultracold neutron losses together with surface wave scattering loss and ultracold spectral evolution, 
a correction to a previous neutron lifetime experiment is found to be insignificantly small. The 
source of the systematic variation between experiments remains to be identified. 

PACS numbers: 
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I. INTRODUCTION 



A recent determination of the neutron lifetime by Serebrov et al. have yielded a value of 
878.5 ±0.7 sta t±0.3 sys s |l|, which is in significant disagreement with the PDG recommended 
value of 885.7 ± 0.8 s [2|. The purpose of this note is to report the results of a theoretical 
study of the effects of surface wave scattering on the gravitational corrections to one of 
the moderate-accuracy experiments [3j that is included in the PDG determination. These 
results might also have implications for the most accurate experiment because the effects 
described here have not been included in the modelling of the experimental system 4j. 

The use of liquid Fomblin, a perfluorinated ether with very low nuclear absorption, to 
produce a nearly perfect UCN storage bottle was suggested and pioneered by C. Bates j^. 
The highest accuracy storage experiments have employed Fomblin coated storage bottles. In 
order to get rid of the thermal scattering, which is a dominant loss near room temperature, 
[l| employed frozen Fomblin and achieved a loss per wall reflection of less than 5 x 10~ 6 , 
a record. However, this is an order of magnitude larger than expected and is a source of 
concern. 

Neutron lifetime experiments employing the storage of ultracold neutrons (UCN) are 
generally based measuring the total storage lifetime, which includes both beta decay and 
wall losses, as a function of volume (or mean free path between collisions) and extrapolating 
to infinite volume. In principle, this seems like a good idea; in practice, however, there 
is a distribution of neutron velocities and each velocity group has a different lifetime. To 
get around this problem, in j^], the volume and storage times were adjusted so that the 
same number of wall collisions occurred for each volume. Because the times were different, 
the neutron lifetime could in principle be directly determined, independent of the specific 
loss mechanism. However, the effects of gravity spoil the simple mean free path scaling 
because the wall collisions are not identical on all the internal storage cell surfaces; given 
that the neutron potential is 1.04 neV/cm in the Earth's gravitational field, the change in 
UCN energy between the top and bottom of a storage cell is about 31 neV, a significant 
fraction of the wall potential, 106.5 neV for the liquid form of Fomblin used in jjj. Recent 
work has shown that surface capillary waves are a significant source of UCN loss due to 



inelastic upscattering beyond the wall potential 



6j. The scattering will result in the change 



of the velocity spectrum of stored UCN. The analysis for the gravity correction that was 
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performed in [3| did not include these effects, for which the gravity correction possibly scales 
differently from the usual loss function. 

Only the minimal published details of the experiment 3| are available. The lack of 
details and inconsistencies between the two publications 3|, |7| is astounding, and rather 
than adjust the model to exactly reproduce the results reported in those publications, a 
different approach was taken. A model analysis was performed with and without the surface 
wave contribution, and the differences between the extrapolated lifetimes for the two cases 
were used to determine new correction terms to be added to those already included in j^j]. 

In this calculation, the UCN energy will be given in cm, representing the change in energy 
with height of a UCN in the Earth's gravitational field. The zero of energy will be take at 
the bottom of the storage cell. Also, the neutron lifetime will be taken as infinity, so the 
correction to 3] will be represented by any residual non-zero loss rate after the scaling 
reduction is applied to the numerical results. 

In jjjj], the UCN storage cell was a box H = 30 cm high, W = 40 cm wide, with a variable 

length L that was used to change the mean free path: 
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A = H- 1 + W- 1 + L- 1 ' ^ 
The UCN energy had a minimum value set by an Al window (53 cm potential) on the 6.4 
cm diameter UCN guide, with a 20 cm rise to the storage cell. The UCN enter the bottle 
at height about 10 cm from the bottom of the bottle, so the minimum UCN energy is taken 
as 43 cm. There in no window on the separate guide that leads to the detector, so the 
minimum energy UCN detected will by 10 cm. The results were insensitive to the assumed 
UCN energy spectrum, so a simple initial 

n(E, 0) = VE for E > 43 cm (2) 

form for the initial density spectrum delivered to the bottle entrance was assumed. However, 
in [?| it is stated that raising the storage cell by 20 cm resulted in a factor of three change 
in the number of stored UCN, suggesting a more significant high energy spectrum. It is 
likely that the large change in the number of stored UCN with height is the result of an 
experimental artifact that is unknown. However, assuming a spectrum of form n(E) = E 2 
did not significantly change the results of the calculations presented here. 

An electrical analog model of the experimental system shown in Fig. 1 will help in the 
following discussion. The storage cell is represented by a capacitor that is charged (filled 



with UCN) by closing the switch located after R in , and connects the capacitor to a current 
source. A sliding door serves as the switch (valve), and in order for the scaling to work, the 
opening area is scaled (e.g., the valve is partially opened) with the storage volume so that 
the ratio of the impedance presented by the filling hole and guide (Ri n ) to the impedance 
representing the storage cell loss (Rb) is constant. In the absence of gravity, 

Rin(E) = R B {E) = 4 (3) 

v{E)A in n tat {E)v{E)A B 

where v(E) = y/2gE, (g = 980cm/s 2 ) is the UCN velocity for energy E, A in is the input 
valve opening area, (fj, tot (E)) is the appropriately gravity average total UCN of energy E 
loss probability per reflection, and Ab is the bottle area (Vb is defined here as the bottle 
volume). It is stated in [3| that the fill hole opening area was varied to keep the loading of 
the source constant. Unfortunately, how this adjustment was performed is not given; one 
might assume that the sum fill hole and bottle impedances was kept constant; this requires 
that the fill hole be reduced for larger volumes. Although this keeps the spectrum constant 
within the guide, the effect on the neutron spectrum in the bottle itself is enormous. There 



is a further contradiction in {7] where it is stated that the filling time was scaled with A. 
Changing the filling time cannot possibly work, so it is assumed here the filling area was 
varied. This effect is referred to as "loading the source" but it seems more profitable to think 
of the effect in terms of the relative impedances represented by the fill hole conductivity and 
the storage cell losses. It should be noted, however, that there is a gravitational correction 
to Rb that spoils the simple scaling. It would make more sense that the fill hole to the 
bottle area, A^/Ab, and that is assumed in the calculations presented here. 

In operation, the storage cell was filled for five or so net storage cell lifetimes, assuring 
equilibrium, after which the fill switch was opened (fill valve closed). After a the storage 
period, the switch to the detector was closed (valve opened) and the UCN counted. It 



is stated in |3j that the guide lifetime was 15 s, but there is no statement of varying the 
opening area of the valve to the detector with storage cell volume. It is stated in [3] that the 
detector valve opening time was varied so that the number of collisions was constant over the 
experiment cycle. Again, details of this time shift, possibly up to 30 seconds or so, are not 
described in any of the publications. The numerical modelling was unstable regarding the 
form of the UCN detection, and shifts in apparent lifetime of up to 40 seconds were possible. 
To get around this problem in the model presented here, the number of UCN left in the 



bottle at the end of a fixed storage period was determined, and as stated earlier, the change 
in number with and without the surface waves was used to determine a correction to the 
correction presented in [3j so the results presented here are independent of these complicated 
and unknown details, and are expected to be a reasonably small percentage effect compared 
to what would be determined with a full model. 



The average loss function is determined by the formalism presented in 
usual UCN loss due to absorption on the walls is given by 



0,3- 
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First, the 



(4) 



where \i is the loss probability due to nuclear absorption and molecular upscattering per 
surface collision, / = 8.3 x 10~ 6 is the ratio of the imaginary and real parts of the surface 
potential, and V is the wall potential, approximately 104 cm for Fomblin, at 283 K [6] . The 
effective volume of the bottle, in the presence of gravity, is given by 



1(E) = jf \]^~ A(h)dh (5) 

where A is the cross sectional area of the bottle at height h, constant and A{h) = LW for a 
rectangular bottle, and we do not assume all neutrons in the storage cell can reach the roof 
(to allow for spectral evolution during storage). The gravity averaged loss rate for UCN of 
energy E is then 

1 v(E) r min ( E ' H "> E -h 
Tw = 4^(E) I » {E ~ k) S{k) dk (6) 



and 



S(h) = LW(S(h) + 5(h- H)) + 2W + 2L (7) 



where S(x) is the Dirac delta function. 

The surface wave upscattering/downscattering can be treated by a similar formalism. 
From {(J, it can be seen that the probability to scatter a single UCN with a given initial 
energy to a final energy bin of width 5Ef is reasonably well described by 

P{Ei -> Ef)5Ef = Eiiane-^-W + a 2 e-^V\ E ^ E f\)5E f (8) 

where a\ = 1 X 10~ 8 /cm, Pi = .065 cm -1 , a 2 = 1 x 10~ 5 /cm, and fa — 2 cm -1 / 2 for Fomblin 
at 283 K. Noting that the energy differences and bin widths are independent of h, the rate 
of change in spectrum can be calculated at h = using the same formalism as used with 



wall absorption to account for gravity, the gravity averaged rate of scattering from Ei to 
EfSEf is 

IW* = i f^MS ^| ^ E^l S(h) x * rfk . (9, 

The spectral evolution from surface waves is then 



dn{E,t) ~ J ^n(E l ,t)T sw (E i ^E)dE i (10) 

/ T sw (E^E f )+ / T SW (E — ► Ef)dEf 
Jo Jv 



dt 

i 

(11) 



n(E,t) 
dn in (E) 



n(E, t) [r SWjOU t r sw> ioss ] (12) 



(13) 

where the first term on the right hand side represents UCN scattered into final energy E 
thus contributing to n(E, t), the second term UCN scattered out of energy E, and the third 
term represents UCN scattered from E to energies higher than V and are thus lost from the 
system. These equations conserve the total number of UCN in the system (total remaining 
stored plus those lost). 

The equations were numerically integrated to determine the UCN spectrum and number 
density as a function of time. The initial spectrum was taken as 

1 in J- sw,loss 

where 

r "<*> = R-kv B < 15 > 
(the spectrum was not assumed to evolve due to the in, out scattering during filling). The 
spectrum was divided into bin 1 cm wide, and integrated in time over 1 s intervals: 

n(E,t + St) = n(E, t) + dnm ^,t) ^ _ ^ ^ ^ + ^ ^ + ^ 
and the total number of stored UCN that can be detected as a function of time is 

N(t) = [ ~f(E)n(E,t)dE. (17) 

Results for the spectral evolution for initially monochromatic UCN of 70 cm energy, with 
and without gravity, are shown in Fig. 2. Results of scaled extrapolation of the neutron 
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lifetime however indicated an effect of less than 0.5 sec on the change in the neutron lifetime. 
Figure 3 shows typical data that was used in this extrapolation. Numerous tests for effects 
of initial spectrum, bottle emptying procedure, and source loading, were performed. It 
might be surprising that such a large effect on the UCN spectrum does not lead to a large 
correction, but insofar as the gravity correction follows the scaling law as used in [3j, the 
effect should not contribute. Indeed, shown in Fig. 4 are the results for storage at time/mfp 
scaled calculations. It can be seen that the effects of gravity are negligible, in that there 
is a less than 1% difference in the spectra, and the majority of the spectrum remains with 
±10 cm of the initial energy of 70 cm. Reasonable agreement with the corrections given in 
{3]] was obtained. In addition, the effects of the surface wave was determined by adjusting / 
for the cases with and without surface wave to give approximately the same loss rate. 

With that said, it was not possible to fully model the experiment described in [3j with the 
available information. This study would be better done in regard to the experiment described 
in 4j , which appears as more complicated, so the simple electrical models used here are not 
applicable. However, there are enough experimental details that exist in Ph.D. dissertations 
andpublished papers to allow construction of a reasonable model of the apparatus described 
in [4|. 
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FIG. 1: Electrical circuit analog of the neutron storage experiment. 
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FIG. 2: Evolution of an initially monochromatic UCN spectrum (70 cm) due to surface wave 
scattering, L = 85 cm. Loss due to upscattering and wall absorption are included, but the beta- 
decay rate is set to zero (as it is for all calculations in this note). 
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FIG. 3: Lifetime/mean free path scaled numerical results. Extrapolation to infinite mean free 
path determines a gravitational correction factor. Effects of finite numerical precision were made 
negligible by considering cases with and without gravity, and with and without surface waves. 
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FIG. 4: Spectrum for two cases: L = 85.71 cm, t = 1000 s compared with L = 12.24 cm (one-half 
mfp), t = 500 s. The differences are small, and can be used to explain the fact that the surface 
wave scattering does not affect the scaling extrapolated neutron lifetime. 
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